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Abstract 



D , We explore spherically symmetric stationary solutions, generated by "stars" with reg- 

, 

ular interiors, in purely massive gravity. We reexamine the claim that the resummation 
^ ■ of non-linear effects can cure, in a domain near the source, the discontinuity exhibited by 
' the linearized theory as the mass m of the graviton tends to zero. First, we find analytical 
difficulties with this claim, which appears not to be robust under slight changes in the 
form of the mass term. Second, by numerically exploring the inward continuation of the 
class of asymptotically fiat solutions, we find that, when m is "small", they all end up in 
a singularity at a finite radius, well outside the source, instead of joining some conjectured 
"continuous" solution near the source. We reopen, however, the possibility of reconciling 
massive gravity with phenomenology by exhibiting a special class of solutions, with "sponta- 
neous symmetry breaking" features, which are close, near the source, to general relativistic 
solutions and asymptote, for large radii, a de Sitter solution of curvature ~ m^. 
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1 Introduction 



Over the past few years there has been considerable discussion about theories with hght 
massive gravitons in their spectra. These kinds of theories first arose in the context of 
brane-world models [1-8] and at the linearized approximation they predicted modifications 
of Newton's constant or even of Newton's law itself at cosmological scales. The evidence 
for a dark energy component in our universe and the associated cosmic acceleration [9, 10] 
made the study of these models rather topical. 

Indeed, generalizing these theories at the non-linear level [11], revealed that they natu- 
rally give rise to a period of late time acceleration of the universe [12-14]. The non-linear 
analogue of a collection of light massive gravitons is a theory with many interacting met- 
rics, one of which only couples to the matter fields of our universe [11]. Depending on 
the particular form of the coupling of the several metrics, the resulting acceleration could 
have interesting testable differences from a (scalar field) quintessence model, as for example 
anisotropic features [14]. 

Despite the interest in these theories regarding cosmology, there are potentially danger- 
ous issues associated with the presence of extra polarization states of the massive gravi- 
tons [15-17]. In particular, the massive gravitons have a scalar-like polarization state whose 
coupling to matter does not depend on the mass of the graviton. This scalar coupling to 
matter is formally analogous (in the linearized approximation) to a Jordan-Fierz-Brans- 
Dicke coupling uj = 0. This coupling modifies (by a factor {2uj + 3)/(2c<j + 4) = 3/4) 
the usual general relativistic relation between interaction of matter and light. If General 
Relativity (GR) is modified by a mass term, or if GR is augmented by the addition of an 
extra massive graviton which dominates^ the matter couplings, the discrepancy between 
the theory of massive gravitons and GR in e.g. the bending of light by the Sun, would be 
at the 25% level, when the current observational accuracy is better that one part in 10^. 

In addition to this (very serious) experimental difficulty, it has also been shown that 

massive gravity has serious theoretical defects [17-19]. As soon as one goes beyond the 

linearized level, massive gravity has six degrees of freedom, instead of the expected five 

(2s -I- 1, for a massive spin s = 2). This sixth degree of freedom is problematic both because 

it represents a jump, with respect to the linearized theory, in the number of degrees of 

freedom, and because its energy has no lower bound. These arguments seem to exclude the 

^In the case where the massive graviton coraponent were sufficiently subdominant there would be no 
observational discrepancy, but also no phenomenological interest for considering such a massive component. 
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existence of dominantly contributing massive gravitons, however small their mass might 
be. It seems, in that case, that gravity is a unique example of an "isolated" theory whose 
massive deformation can be excluded to infinite accuracy. 

One could question several points in the arguments of [15-17] recalled above. The first 
one is that these arguments are based on the study of models containing an explicit breaking 
of the (linearized) gauge invariance 6h^i, = d^^^ + dy^^ by a Pauli-Fierz mass term for the 
light graviton(s), which reads: 

2 

- — {h,yh^^-e) (1) 

One could hope that the peculiar result of the discontinuity is linked to this explicit break- 
ing and that continuity in predictions might be restored if the mass for the spin-2 field were 
generated spontaneously. For example, if light gravitons arc generated by the compactifi- 
cation of a higher dimensional theory, one could imagine a higher dimensional mechanism 
which would restore the continuity in the theory^. We will not explore such possibilities in 
this paper, but study instead theories with explicit mass terms with the aim of clarifying 
whether such theories are indeed physically sick. 

A second potentially weak point of the above argument is that the effects of Pauli- 
Fierz mass terms were considered only around flat space backgrounds. One could imagine 
that the discontinuity that was found is a peculiarity of just the fiat background and that 
if one considered a background with curvature, some of the difficulties might be evaded. 
Indeed, it was found that for instance in constant curvature backgrounds [{A)dS spaces], 
the extra polarizations of the massive gravitons have a coupling ^ m/H where m is the 
mass of the graviton and H the "Hubble" constant of the {A)dS space [21-24]. In that 
case, the predictions of the massive theory were indistinguishable^ from the massless one as 
long as m <^ H. This, however, prevents the massive gravitons from being cosmologically 
interesting since their Compton wavelength is much bigger than the "Hubble" radius of the 
"bare" cosmological constant. 

In constant curvature backgrounds the parameter controlling the smoothness of the 
hmit is m'^/TZ where TZ ex H"^ is the curvature of the (A)dS space. One might conjecture 
that the same is true for any curved space and that the smoothness is controlled by a 
parameter ~ m^/TZ, where TZ is some curvature invariant. In the case of generic vacuum 
^For an alternative way of generating spontaneously masses for gravitons see [20]. 

^Notc that at the quantum level the discontinuity formally reappears [25,26], but as a quantum effect 
is suppressed and unobservable in any conceivable near future experiment. 
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Einstein backgrounds (with vanishing scalar curvature and Ricci tensor) TZ might be defined 
as some invariant of the Riemann tensor (such as the square root of the square of the 
Riemann tensor). In the case of a Schwarzschild background, TZ ~ Rs/R^, where Rs is 
the Schwarzschild radius of the source. Then a smooth limit might be expected to exist 
for distances R <^ Rp, where the upper bound Rp of the expected interval of existence of 
smooth solutions is: 

= {mRlr/' ^ ^^-^^^''^ 

where = Tn^^ denotes the Compton wavelength associated to the mass m. 

A similar conjecture was made a long time ago by Vainshtein [27] based, however, on 
a different argument. That paper made two basic points: (i) it questioned the use of 
perturbation theory (and especially of the linearized approximation) in the derivation of 
observable consequences of massive gravity by showing that non-linear effects, proportional 
to a negative power of m^, are important in a wide domain around the source, and (ii) 
it conjectured for the first time that the resummation of non-linear effects might actually 
restore continuity in a domain near the source. More precisely, this paper sketched the 
construction of solutions which can be expressed, at least within some intermediate range 
of distances Rg <^ Rq < R <^ Ry, as a series in positive powers of and whose leading 
term coincides with the one of GR. The interval where this expansion was constructed was 
bounded on the left by the radius of the star, it!©, and on the right by the length scale: 

= (^^^^ ^ (A^^5)^/^ (3) 

which differs from the one obtained by the conjecture made in the previous paragraph. 

In the simple example of a massive graviton with Compton wavelength of the order 
of the universe's horizon and of a source as massive as the Sun, the distance scale (3) is 
much bigger than the distance scales on which relativistic predictions of GR are tested. 
Thus, if the conjecture were correct, the massive graviton proposal would not be excluded. 
[Actually, the same holds if continuity is restored on the scale (2), though Rp <^ Ry-] The 
arguments of [27] were repeated, with more intermediate details, in the recent paper of [28] . 
This paper also provided a cosmological analogue of the continuity-discontinuity interplay 
in the context of the branc-induced gravity models (see also [29] for a cosmic string example 
for the same kind of brane models). 
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However, the arguments of [27, 28] are questionable, and one of the main thrusts of 
the present paper will be to show that they actually are inconclusive. The questionable 
aspects concern both the local existence of continuous solutions, their global relevance, and 
even their global existence. For instance, we shall show that the action principle postulated 
in [28] does not admit , in the intermediate range Rs <S Rq < R <^ Ry, local solutions 
of the type conjectured in [27] for restoring continuity. Moreover, it is not at all clear if 
such foca^ solutions, when they do exist (for some other action principles), can be globally 
continued to the asymptotically fiat ones of the linearized theory, a concern first raised 
in [17] and confirmed later in [30]. Actually, it is not even clear if these solutions are indeed 
solutions, i. e. if they do not contain naked singularities at some finite distance. 

Let us note here that there have been several papers [31-35] where the Vainshtein 
idea was investigated in brane-induced gravity for static spherically symmetric sources and 
novel phenomenology was proposed [34, 35] ''. In our view, the approximations made in 
these papers, although plausible, do not address the global properties of the solutions. We 
will not explore the specific brane-induced gravity model in this paper. 

In the present paper we will discuss, within the context of purely massive gravity, the 
claims of [27] as well as the more expanded version of them that appeared in [28]. We 
will analyse in detail the procedure that was used in [27] to obtain the results, shedding 
light to some obscure points, and repeating the same calculation for several mass terms. 
We will show that the method has serious limitations and that there are contradictions 
between the statements made in [27] and in [28]. Indeed, we shall point out that the 

field equations studied in [27] do not derive from the action principle postulated in [28]. 

ActTially. they cannot be derived from any action principle. We shall show that, if one 

starts from the simple action principle written in [28], there exist no expansions of the 

type postulated in [27] as a way to restore continuity. However, we shall show that such 

expansions do exist for other action principles (and do exist for the field equations written 

down in [27] , after correction of some misprints) . We will also see that the scale Ry up to 

which the expansion of [27] makes sense, is not universal for different potential terms. We 

will then study (numerically) the global aspects of such solutions. Our numerical simulations 

strongly suggest that: (i) the Vainshtein-type local (approximate) solutions do not match, as 

they were supposed to, the asymptotically fiat approximate solutions generated by normal 

^See [36] for a demonstration of the discontinuity at the linearized level for spherically symmetric sources 
in brane-induced gravity models. For the same models, see also [37, 38] for special spherically symmetric 
solutions. 
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perturbation theory (in agreement with [30]), and in fact that (ii) the asymptotically flat 
solutions cannot be extended inwards to globally regular solutions. Even if we forget about 
the issue of mass-continuity, we find the striking result that all asymptotically flat solutions 
run into naked singularities as the radial coordinate R decreases. 

However, our conclusions will not all be negative. Indeed, we shall explicitly construct 
some globally regular solutions of massive gravity which are continuously connected, when 
m — > 0, to GR solutions, and which are phenomenologically consistent with experimental 
tests of relativistic gravity. The global solutions we shall (numerically) construct are gen- 
eralizations (to the inclusion of a central source) of the black-hole-type solutions of massive 
gravity constructed long ago by Salam and Strathdee [39]. [The latter solutions were gen- 
eralized by Isham and Storey [40] for a particular class of bigravity theories [41].] The main 
difference between these solutions and the ones studied in Refs. [15-17,27] is the behaviour 
at infinity. The solutions are not required to be asymptotically fiat, but instead to match to 
a cosmological solution of massive gravity (which is de Sitter, in the case at hand). In fact, 
a general argument of [11] (see section 4 there) has shown how to construct (at least for a 
limited time) general classes of solutions which arc continuously connected, when m ^ 0, to 
GR solutions representing local gravitating systems (such as the solar system), embedded 
in some global cosmological background. Another big difference with the solutions stud- 
ied in [15-17,27], and in the first part of this paper, is that these cosmologically-matched 
solutions are of a special "symmetry- breaking" type (see below) . 

As we shall discuss in our conclusions, our results leave open several important issues 
which must be tackled before a firm conclusion can be reached concerning the physical 
consistency (or inconsistency) of massive gravity theories. 

2 Action and potentials for massive gravity 

Our starting point is a generic action for (four dimensional) massive gravity. This can be 
obtained, as it was shown in [11,39], from a four dimensional bigravity action if wc send 
the gravitational constant associated with the second metric (the one not coupled to "our 
world") to zero. Then the corresponding metric gets (formally) frozen, i.e. it becomes a 
non-dynamical Einstein space background. We assume that this non-dynamical metric is 
flat, and we accordingly denote it by f . [We do not necessarily assume that the metric f is 
written in Lorentzian coordinates; i.e. the components f^iu{x) are not necessarily assumed 
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to be simply 77^,^.] The only remaining dynamical metric is denoted g. Then the effective 
action for the purely massive gravity is: 



-gR[s] 



m 



(4) 



where G is the "massive" version of Newton's constant and m a mass parameter, which 
for the specific normalization of the potentials V which we will consider, is the canonically 
normalized Pauli-Fierz mass of the graviton. In the following, we explore the following 
possible potential terms: 



^{tr[((g-f)OV(tr[(g-f)nf} 

7{tr[((g-f)f-^)V(tr[(g-f)r^])^} 
^{tr[((g-^ - Of)V (tr[(g-^ - f-^)f])'} 
7{tr[((g-^ - r^)f) V (tr[(g-^ - Of])'} 



(5) 
(6) 
(7) 
(8) 



Note that the above quantities are not scalars, but scalar densities. Comparing with the 
notation of [11], V = (density) F where F is a scalar. All of the four above potentials are 
of the Pauli-Fierz type near g f . Note that the first one was written to be the starting 
point in [28]. 

Each contribution to the action gives rise to a corresponding contribution to the 
energy-momentum tensor: T^^ = —(2/-^/^) SSx/Sg^". The "gravitational" energy-momentum 

tensors arising from the mass terms will be denoted Tjf} . For the specific mass terms above, 
they are respectively: 





m 


16nG 






IGTrG 




2 r 

m 


IGTrG 




2 

m V 


IGTrG X 



-9 



{9fiK9u\h'^ -hg^^^g^xf 



-f 



(9) 
(10) 

(11) 
(12) 



where we have defined h^i, — g^^, — f^^, and H^^ — g^" — f^^. The indices of h^^, are raised 
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by f^"" and the ones of H^"' are lowered by /^^. Note that H^''' = -h^"" + h'^'h''^ + ■■-, which 
explains the sign differences in the above expressions. 

In this paper wc wish to consider spherically symmetric stationary (SSS) solutions for 
these mass terms, i.e. SSS solutions of: 

R,. - \r9,. = SnGiT^i^ + T^-^")) (13) 

Note that the separate diffeomorphism invariance of the matter action 5matt (which is 
assumed to couple only to g) implies the separate conservation (on matter shell) of the 
material energy-momentum tensor: 

A consequence of Eqs. (13), (14) and of the Bianchi identities is then the separate conser- 
vation of the gravitational energy tensor: 

V^T^^) = (15) 

3 Poincare-covariant perturbation theory 

Note that all the actions (4) admit as exact solution, in the absence of matter, the "trivial 
vacuum" g = f. [See, however, section 8 below for a discussion of non-trivial vacua.] 
When representing the non-dynamical flat background metric f in Lorentzian coordinates, 
^■e. fnu — Vtiu, one can then develop a Poincare-covariant perturbation theory: /i^^ = 

Qfiu — = + h^^} + ■ ■ ■. Within such an approach it is implicitly required that the 
massive gravitational field h^i, decay at infinity, so that: 

g^lu ^r]^^ , as r ^ oo (16) 

In the first part of this paper, we shall impose this standard requirement, the problem being 
to assess the existence and continuity of solutions of massive gravity matching the trivial 
vacuum at infinity. [However, we shall relax this requirement in the second part of the 
paper.] 

Before entering the details of our investigation of non-hnear SSS solutions, let us, as a 
warmup, recall the basic features of perturbation theory. At the linearized approximation. 
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the equations of motion read: 



Here all index raisings and lowerings are made with 7]^!^, and the source term T^j, is equal, 
at the linearized approximation, to the matter energy-momentum tensor. [See below, for 
the generalization of Eq. (17) to the non-trivial vacuum cases.] Following [17] we have 
introduced a parameter a to study the special role of the Pauli-Fierz mass term, which 
corresponds to ctpp = 1. 

Let us recall the consequences of Eq. 17). The divergence of this equation yields, when 
using the conservation of the source, d^T^,, — 0, the constraint: 

Then the trace gives: 

2(l-a)n/i + m^(l -4a)/i= 167rGT (19) 

Prom the latter equation we see that if a = 1, i.e. when we sit on the Pauli-Fierz point, 
the trace of hf^i, is locally determined by a constraint: 

and thus substituting in Eq. (17) one gets: 

1 , 1 /„ 1„ \ 1 1 



This result for the massive gravity field in terms of the source exhibits the difficulties of 
massive gravity. First, there are the phenomenological difficulties associated with the 1/3 
instead of the Einsteinian 1/2 factor in front of Trj^j^^. Second, there are the theoretical 
difficulties associated with the presence of a factor mT^ in some terms of the solution. 

Before discussing further these difficulties beyond the linearized level, let us also briefly 
recall the results of [17, 19] concerning the non-Pauli-Fierz mass terms, i.e. the case a ^1. 
In this case, we see from Eq. (19) that the trace of /i^,^ becomes a dynamical field, i.e. a 
new degree of freedom. Then the corresponding result for the expression of h^^ in terms of 
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the source is: 



1 



1 



— □ + 




2a -1 1 1 




(22) 



2(1 -a) -□ + m2-n + m; 



where 




2 4(3; - 1 
' 2(1 -a) 



(23) 



is the mass of the extra scalar degree of freedom. [It is non-tachyonic if a is in the interval 
1/4 < a < 1]. Note that this field is always ghostlike (independently of the value of a), 
as seen by considering the various contributions to the action, /i^j.T^'^, where the term 
— {1/6)1] niy{—n + mQ)^^T contributes with the opposite sign from the first term (note that 
the last, double gradient term, does not contribute, after integration, because of the conser- 
vation of the source, d'^T^^ = 0). As emphasized in [17, 19] the case a ^ 1 does not exhibit 
the difficulties of the Pauli-Fierz case. On the one hand, the expression (22) contains no 
dangerous denominators, vanishing with m^, and on the other hand, it smoothly merges 
(thanks to the identity 1/3 + 1/6 = 1/2) into the GR result for the integral of hniyTi^" when 
— > 0. Those nice continuity properties of the case a 7^ 1 pointed out in [17, 19] are 
particularly evident in the case a — 1/2, but hold in all cases a ^ 1 , independently of 
whether the extra scalar degree of freedom is tachyonic or not. [These continuity properties 
have been recently further studied in [42,43]]. Nevertheless, the presence of a new ghosthke 
degree of freedom makes the theory for a 7^ 1 pathological as a quantum theory. In view 
of this quantum ghost instability, we restrict ourselves to the ghost-free Pauli-Fierz case. 
We note in passing that Pauli-Fierz mass terms naturally arise from higher-dimensional 
gravity models. Indeed, the structure of the Pauli-Fierz mass term is already encoded in 
the structure of the gradient terms in Einstein's action (see, e.g., section 3.1 of [11] ). 

Let us now briefly indicate how the Poincare-covariant perturbation theory would pro- 
ceed beyond the linear approximation. As is done in the Poincare-covariant perturbation 
theory in GR one can define an effective energy tensor T^f, which combines the matter 

energy tensor with the non-linear parts ~ ddhh + irP'hh + ■ • ■ of the left-hand side of the 
exact field equations (13), such that the full field equations read as the linearized ones, Eq. 
(17), with the replacement T^^, —>■ T^. This effective energy tensor is not exactly conserved 
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but rather satisfies a relation of the type d^T^ — m^{dhh H This modifies both Eq. 

(18) and Eq. (19), and thereby Eq. (21), by non-hnear terms. By iteration, one can then 
deduce the following structure for the non-linearity expansion of the solution: 

Here, we have denoted by C/ a typical gravitational potential which is, outside the source 
and for distances i? <^ A^, of order U Rg/R. The expansion (24) contains increasingly 
high inverse powers of m. We can rewrite it as a sum of terms of the type {Rp^q/ R)"' with 
a sequence of length scales defined as: 

it:,,, ^ (A^^i?|)V(p+^) (25) 

where, as above, = m~^. 

The expansion above can numerically make sense only if the radius R is (much) larger 
than all the relevant length scales Rp^q. In view of the extremely large difference between 
A^ (~ 10^^ cm, say) and Rs (~ 10^ cm, say, for the Sun) the various length scales tend 
to be extremely large, which confirms the first point made in [27], namely the irrelevance 
of perturbation theory for describing the massive gravitational field near the source. For 
instance, the light deflection by the Sun probes the gravitational field just outside the radius 
of the Sun, while the expansion above represents h there as a sum of terms, starting with a 
small first term U ~ Rs/R ~ 10~^, but continuing with extremely large "corrections" such 
as V/imRf ~ 10+22 I 

Actually, the straightforward perturbation expansion (24) can be significantly improved 
by gauging away some of the worse terms. Indeed, the primitive source for the "bad" inverse 
powers of is the last term in Eq. (21). However, this term is mainly a "gauge term", in 
the sense that a suitable coordinate transformation x^j^ x^g„ can remove this term from 
the physical "massive" metric g^uix), at the cost of introducing it in the non-dynamical fiat 
metric f^uix) (now written in the "curved" coordinate system x'^^^ instead of the original 
Lorentzian coordinate system). However, the crucial point is that this "gauging away" 
cannot remove all inverse powers of rri^. Indeed, if we write x'^^^ — x'^^^ + C^{xnew), we 
have g^^{x) — dij,{x°' -\- ^°')d,^{x^ + i^)ga^{x -\- ^), so that the expansion of the new field 
^new = g,new _ ^ contains, bcsidcs the terms which can remove the last term in Eq. (21), 
the new terms ~ d^h -\- d^d^. Therefore, with ^ ~ dU jvr? ~ U/{rin?R) chosen to remove 
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the worst term oc m ^ in Eq. (21) , we still end up with several terms that contain inverse 
powers of m^, namely: 

One could continue to apply suitable gauge transformations to remove some of the worse 
terms arising from computing higher iterations of h, but it is clear that once some negative 
powers of have entered h they will keep generating higher powers in higher iterations. 
In view of the vast difference in scale between ( ~ 10^^ cm) and phenomenologically 
tested scales R (~ 10-*^^ cm for the Sun radius) the dangerous inverse denominators are of 
order of powers of {rriR)^^ ~ 10^"^^ and completely outgrow the small numerators which 
arc powers of f/ ~ 10~^. Finally, the improved expansion (26) is again irrelevant for 
describing the massive gravitational field near the source, and can only he valid in some 
neighbourhood of infinity. On the other hand, note that if we consider, for simplicity, the 
case of stationary sources the basic "Newtonian" potential U in the expansions above will 
always be a Yukawa-type potential (obtained by the action of {—\3 + m^)^^ on the material 
source T). This potential will decrease exponentially as i? ^ oo. Therefore it is clear 
that perturbation theory will be formally well defined, to all orders, when R oo, and 
will generate a numerically meaningful scries (say in the whole domain R > m^^). This 
series can be thought of as uniquely defining a specific solution of massive gravity. The 
main question of concern being to know whether this solution, initially defined only in a 
neighbourhood of infinity, can be continued into a regular solution everywhere. 

Before leaving the topic of Poincare-covariant perturbation theory, it is important to 
notice the following. If we apply this perturbation theory to the particular case of static 
spherically symmetric sources (with no motion in the Lorentzian coordinate system), it is 
clear that each term of perturbation theory will be static and spherically symmetric. In 
particular, each term will be invariant under time-reversal t —t. This implies that, in 
Lorentzian coordinates (i.e. with f^^, — 1]^,^,) there cannot be off- diagonal terms goi in the 
massive metric. Therefore, the exact solution defined (as above) by perturbation theory 
can be written in a bi-diagonal form, say: 

i = -dt^ + dr^ + r^dn^ , g = -a^rft^ + ft^dr^ + cW, (27) 

with a, 6, c some functions of r, and dfl"^ — dO"^ -\- s\v? Odcf)'^. We shall later reconsider this 
bi-diagonality. 
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4 The A, V gauge 



When discussing SSS solutions, it is instructive to work in two different gauges for the 
metrics. We have just introduced the a, &, c gauge (27) (which will be used below to intro- 
duce some convenient variables c, c, 6). It is also convenient to work in the A, /x, v gauge 
introduced now. The latter is useful for comparing our results with the literature, while 
the former is similar to the gauges used in cosmological studies, and has also the merit to 
give a convenient Lagrangian formulation to the system. Let us write the observable line 
element ds^ — g^^dx'^dx^ and the reference flat metric one ds^ = f^i.dx'^dx'^ in the following 
( "Schwarzschild" ) gauge: 

= -e'^^^Uf + e^^^'UR' + R^dn'' (28) 
dsl = -df + (^1 - ' e-^^^'^dR' + e-^^^'^R'dn^ (29) 

where ' = This gauge was used in part of the literature of spherically symmetric 
solutions in massive gravity. It has the advantage of separating the directly observable 
gravitational variables u{R), X{R) from the "gauge" function //(-R). The coordinate redef- 
inition r = Re~^^'^ shows trivially that the second metric is indeed flat, and that the form 
(28), (29) is actually equivalent to the form (27). [See below for the explicit hnk between 
the two sets of variables.] 

The most general matter energy-momentum tensor which respects the spherical sym- 
metry is of the form: 

Vtt) = diag(-p, P„ P„ P,) (30) 

where we have left open the possibility of having different different radial P^ and tangential 
Pt pressures. In the case of a "fluid" source, we shall have P^ = P*. Then the (t,t) and 
(P, R) components of the Einstein equations read: 

= m^ft{X,ii,u,iJ,',R) + 8TiGpe'' (31) 

^ + ;^(1 - e^) - m^fRiX, fi, i^, fi', R) + SttGP, (32) 
where the quantities ft and fu are proportional to the {t, t) and (P, P) components of the 
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energy momentum tensor T^fJ generated by the potential term, and are defined as: 



1 1 5V _ SttG (,) 1 1 5V _ SttG (,) 

From the conservation of the energy-momentum tensor generated by the mass term, 

i.e. from VT^k^^O, we obtain a third equation from the K — r component (the k — t,9,(f) 
components are identically zero). Defining: 

/» - ^V"T(;) (34) 

lit 

we have the following constraint equation for m 7^ 0: 

fg{X,fi,i^,X',fi',i^',fi",R) = (35) 

Note that this equation has no m dependence. In the Appendix A we present the functions 
ft, fn and fg for all the potentials that we are considering. 

Finally, from the conservation of the matter energy-momentum tensor we have the final 
equation which closes the system: 

2P; + 1(3P, - p - 4Pt) + (^u' + ^{P + Pr) = (36) 
Note that in the "fluid" case, i. e. for an isotropic pressure, the latter equation simplifies 



to: 



P'-~{P + P) (37) 



4.1 The General Relativity m = limit 

In the GR limit, one must discard Eq. (35), and set to zero in Einstein's equations 
(31), (32). It will be helpful to remind ourselves of the solution of the Einstein equations in 
GR for a star of constant density p = po and isotropic pressure P = Pr = Pt- In this case 
(considered here for simplicity), the exterior (vacuum) solution , after absorbing a constant 
by a time rescaling, reads: 

e^^^l-^ , X^-u (38) 
13 



where Rs = 2GM is the Schwarzschild radius and appears here as an integration constant. 
On the other hand, inside the star, eq. (31) can be integrated to: 



e-' = l- ^-^^ (39) 

where we have defined m{R) = ^irpoR^. Note that we have imposed the constraint \{R = 
0) = 0, which is necessary to avoid a conical singularity at i? — > 0. If wc denote the 
radius of the star by Rq, matching to the external solution gives the links m^Ro) — M and 
Rs — ^ttGpqRq. We then can solve the remaining equations and obtain the well known 
result for the pressure: 




and the redshift function u (corresponding to the same rescaled time as the external solution 
above): 




4.2 The perturbative limit for m 7^ 

Let us now see how the results change if we consider the massive theory. We will firstly 
deal with the linearized equations of motion with the aim of understanding analytically 
how one can match the unique decaying exterior solution into a regular interior one. In 
other words we will assume that X,iJ>,iy <S 1. We will solve these linearized equations in 
the interior and the exterior of the star and then match them to determine completely the 
various integration constants. 

4.2.1 The exterior star solution 

The equations (31), (32), (35) at linear order in A, are the following (for all four mass 
terms) : 

A' A /X / X 

^ + ^--^i>^ + ^l^ + Rl^) (42) 
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(43) 
(44) 



Note that the hnearization of the equations has removed entirely the function from the 
constraint equation (35). These equations, give a simple second-order differential equation 
for V. 

v" + '^v' - = (45) 
R 

Then the other variables are obtained from v via: 

R , v' V , , 

The generic solution of the v equation above contains two integration constants, C\ and 
C{. The solution proportional to C\ is exponentially decaying at infinity, while the one 
proportional to is exponentially growing. In keeping with our boundary condition (16) 
we shall reject the exponentially growing solution and consider only the other one, namely: 

V = -^e-™^ (47) 
R 

A=^fl + ^V-"^^ (48) 

2 \ mRj ^ ^ 

Note that, if needed, the exponentially growing solution is simply obtained from the lat- 
ter solution by changing m —m. The integration constant Ci is proportional to the 
Schwarzschild radius Rs = 2GM of the source (the precise link will be seen later when we 
match with the interior). In the region R ^ we obtain the limits of the above solution 
for A, II, v. 

Ci V Ci Ci , , 

Note that the second relation differs from the corresponding GR relation A ~ —u. These 
show the well known fact that there is a finite discontinuity between this metric and the one 
obtained in GR, independent of the mass of the graviton as was shown in [17]. Additionally, 
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the last relation shows that, in the X, iJi,i^ gauge, and at the linearized approximation, there 
is a dangerous denominator only in the "gauge" function /i. Actually, one can see that 
the gauge function /i is similar to the gradient of the coordinate transformation that we 
introduced above to remove the leading 0{m~^) terms in the (linearized) perturbative field 
h. We had above ~ U/{mR)'^ which indeed corresponds to ~ Rs/{m'^R^). Therefore, 
as above, the choice of gauge has allowed us to "regularize" the linearized physical metric 
u, A, but we expect that non- linear effects will re- introduce some negative powers of in 



4.2.2 The interior star solution 

In order to be able to apply our perturbative calculation to the problem of source matching, 
we need to formally consider a star with sufficiently small density, and sufficiently large 
radius, so that our assumptions A, ^ 1 remain satisfied everywhere. We shall give later 
(after having taken into account non-linear effects) the precise conditions that the star 
characteristics must satisfy. Evidently, such a low density "star" is not physically relevant 
to the discussion of real stars, but the point of this subsection is to show in detail how, in 
principle, one can match the well-defined exterior perturbative solution to a unique interior 
one. 

The interior linearized equations of motion (again for all four mass terms) are: 



u, A. 



R^ R^ 




(51) 




(52) 




(53) 



and additionally the continuity equation is: 




(64) 



These equations, give a differential equation for u: 




(55) 
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Let us assume, for the sake of simplicity, that p — const. = po and P <^ po, i.e. a constant 
density star of low pressure, so that we can neglect the pressure P in the above equation. 
This assumption means that we are considering a non-compact start, with Rq ^ i?5, i.e. 
Rq <^ (SttGpo)"^^^- We can then safely neglect the second addendum in the u' parentheses. 
Thus, we have the simple differential equation: 

u" + - m^v = SttG^Po (56) 
R 3 

which differs from the exterior equation (45) above by the source term on the right-hand 
side (RHS). Under the same approximations, the other metric functions are determined in 
terms of a solution of the interior v equation by means of: 

A = , [X = - - (57) 

The general solution for the v function that is regular at i? = depends on only one 
integration constant, say Ci , and reads: 

.(i?) = -^^ + ^sinh(mi?) (58) 
6 R 

The continuity of v relates the two integration constants, C\ of the exterior solution and 
C2 of the interior one, as: 



Rq 



sinh( mi?©) 



4 SttGpo _ C^^-mRQ 

3 Rq 



(59) 



The continuity of the first derivative of i/, determines the constant C\ in terms of the mass 
and the radius of the star: 



4i?s 

C\ — - — ——^[mi?0 cosh (mi?o) — sinh(mi?o)] (60) 



Consistently with our problem of studying the m — > limit, we shall assume that the radius 
of the star is such that Rq <^ wT^ . We then get that: 

Ci « ^Rs (61) 
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[In the opposite limit, m -"^ <^ Rq we would have obtained that Ci ft; (mfi|)2 e'"'^°]- Note 

that the factor | here is linked to the well-known fact that, if one assumes the validity of 
linearized theory, one must renormalize the "bare Newton constant" G appearing in the 
massive gravity action by |G = Gn to recover the usual Newtonian constant G^. 

Finally, the full solution is obtained by inserting the matched solution for i/ into the 
expressions (57) above for the functions A and /i. 

Let us finally sketch the structure of the (unique) solution of perturbation theory, in 
the A, II, V gauge, when taking into account the next order in perturbation theory. If we 
introduce U ~ e~"^^Rs/R ~ Rs/R (in the region Rq < R <^ "t."^), we can write: 

u^X^U+(l + ^ + --- + ^^)u' + --- , (62) 
\ niR [mRf J 

^ ~ + ir^ + • • • + + • • • ' (63) 

[mRy \{mRy [mRy J 

This expansion is related to the general perturbation expansions (24), (26) above. In 
fact, roughly speaking the expansion for the "physical variables" A corresponds to the 
gauge- improved expansion (26), while the expansion for the gauge variable corresponds 
to the expansion of the gradient of the "improving" gauge transformation ^ introduced 
in section 2 above. 

One would need to study more carefully the structure of higher order terms in the 
expansions above to delineate what are the most relevant length scales Rp^q , see Eq. (25), 
determining the range of validity of the A, ^i, v perturbation theory, i.e. for determining the 
constraints on Rq and po ensuring that the series above (considered both in the interior 
and the exterior) make sense. As it is clear that the crucial powers of mR entering the 
successive monomials f/"/ {mRy will increase linearly with the perturbation order n, there 
will be a finite limit, when n ^ oo, to the sequence of relevant scales Rp,q, and therefore 
this sequence will have a finite, global least upper bound, or supremum, i?sup- Perturbation 
theory is then valid (everywhere) if Rg^p ^ R&- 

It is evident that since we need (for being able to use perturbation theory) to constrain 
ourselves to stars of extremely small density, the solution that we have obtained is unreal- 
istic. However, we have gone through this explicit derivation to show how, in principle, our 
various physical requirements (decay at infinity, matching at the star radius, and regularity 
at the origin) determine a unique solution in terms of the equation of state (po) and the 
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mass of the star, and the mass of the graviton. Note that this uniqueness of the solution 
is a non-trivial consequence of the treatment based on perturbation theory. Indeed, per- 
turbation theory allowed not only to provide boundary conditions at infinity that selected 
one solution from the second-order differential equation for v, but it also allowed one to 
kill the "sixth degree of freedom" hnked (in the A, u gauge) to the fact that the equation 
fg — is a second-order differential equation in /i. The perturbative treatment allowed us 
to recursively determine the higher derivatives of /i from the algebraic calculation of /i in 
terms of u and u' (see Eq. (46)). 

5 On a claim by Vainshtein 

In the old paper of Vainshtein [27], it was claimed that continuity in the massless limit 
could be obtained if one, instead of treating the problem perturbatively, used a different 
"non-perturbative" expansion, proceeding in positive powers of m^, whose first terms for 
the physical variables A and u were their usual GR values, and whose first term for the 
gauge variable /i were (in the exterior region) of order: 



The aim of this section is to critically reexamine this claim, to explain in detail some points 
which have remained unclear, and to discuss its serious difficulties. 

A first point to clarify concerns the theoretical framework chosen to "define" massive 
gravity. Indeed, the starting point of the original reference [27] was not an action, but 
the Einstein equation with a postulated energy-momentum tensor for a massive graviton. 
However, it is easy to see (by checking the dissymmetry in {x, y) of the functional derivative 

of the postulated g{x)Tj^\x) with respect to g°'^{y)) that the field equations postulated 
in [27] cannot he derived from any action (and, in particular, they correspond to none of 
the four models written down at the beginning of this paper). This makes them quantum 
mechanically inconsistent. The energy-momentum tensor postulated in [27] reads (after 
correcting the sign of the terms as they appear in Eq. (1) there, which actually corre- 
sponded to a tachyonic mass term) : 




(64) 




2 

m 



hr]i^u) 



(65) 



IGttG 
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Note also that, in the same paper, the exphcit (t,t) and (R,R) components of the 

Einstein equations appear with the correct sign, consistent with the T^^^ given above. For 
this energy momentum tensor the functions ft, Jr and fg are given in Appendix A. Note 
also that in [27], fg was written down incorrectly, nevertheless the final result of [27] starting 
from Eq. (65) is mathematically correct. 

The basic idea of [27] was the following. The root of the discontinuous behaviour 
between massive gravity and GR is the "new" constraint fg, whose non-linearity expansion 
starts with (when keeping terms relevant when Eq. (64) holds): 

1 A u' Anil' 7/x'2 ^ ^ , . 

The ellipsis in (66) contain quadratic terms 0{iJiv -\- jiX -\- v'^ -\- ■ ■ ■) as well as higher non- 
linearities. 

When this equation is solved in the usual perturbation theory, i.e. by solving for the 
linear terms 0{\) and 0{v') and then adding the quadratic terms 0{ii^ -\- jiv -\- ■ ■ ■) as 
corrections, one gets the "discontinuous" result (50) above. The observation then was 
that one might cure the discontinuity by assuming that the 0{ii^) terms in Eq. (66) are 
comparable to the linear ones 0{X) and 0{u'). In more detail, one assumes that the three 
functions A, n, v admit an expansion in positive powers of m^, say: 

oo 

f{R) = J2m'yn{R) (67) 

n=0 

where the leading terms, uq, Aq, fio, are such that the first two are solutions of Eqs. (31), 
(32) with set to zero, i.e. (in the exterior): 

-.„W = V«) = -l„(l-f).f.i(f)%... (68) 

while /Iq is a solution of the constraint equation (66) with the above given uq, Aq substituted, 
i.e. 

Ao ly'o 4/io/i[) 7//(2 , . 

^-2^ + ^ + ^ + ^°^° + --- = ° ^''^ 

One is looking for a solution //q of Eq. (69) which satisfies (in the exterior region) Eq. 
(64) at lowest order. Note that the integration constant associated with uo has been set to 
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zero so that one is consistent with the "trivial- vacuum" (see section 3) . Note also that the 
quantity Rs entering Eq. (68) is, at this stage, a GR-like integration constant whose exact 
link with the mass of the source is not important. 

One then finds that the above requirements lead to a unique (exterior) solution for /xq, 
which further admits an expansion in powers of Rs/ R starting as: 




I ^ , , I Rq 39281 R,Q , ,_ , 

A.0=WT^^/^|l + T^Vi + Y3i0ilf + -I (TO) 

where the higher order terms are obtained by considering the constraint equation (69) 
augmented with higher-order terms in Aq, i^q and /iq- 

Starting from this result for /iq (where all higher-order coefficients are, in principle, 
uniquely determined), we then go back to Eq. (31) (in the exterior region), which yields a 
linear first-order differential equation for the C(m^) term Ai in A: 

1 \,, 5 [R^ 257 Rs 



The general solution of this inhomogeneous linear ordinary differential equation (ODE) can 
be written as the sum of a particular inhomogeneous solution: 



™^A. = (™«)M-^A^f-i§ + ...l (72) 

plus the general solution of the homogeneous equation, which is of the form c\i/R. The 
latter general homogeneous solution can be absorbed in a 0{rn?) change of the integration 
constant Rs, and can therefore be physically ignored. 

Wc can then insert the known values of z/qi ^Oj A^o and Ai in Eq. (32) to get a linear 
first-order differential equation for the 0{m^) term i>i in i/: 



R^--T^ = yr3\lT-YAT+--- ^'^^ 

Again, the general solution of this inhomogeneous linear ODE can be written as the sum 
of a particular inhomogeneous solution: 
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plus the general solution of the homogeneous equation, which is simply a constant Ci,i. The 
latter constant can be absorbed in a rescaling of the time variable, and can therefore be 
physically ignored. 

The next step is more tricky. Indeed, we get a second-order hnear ODE for the C(m^) 
term fj,i in m of the type: 



Here a — 7/4 and b — A are the coefficients appearing in the //-quadratic terms in the 
constraint (66). The problem now is that the general homogeneous solution of Eq. (75) 
a priori introduces two new integration constants Cmi+, Cmi- in a solution of the type 
= c^i-^-R^+ + c^i^R^- where the exponents s+, s_ are the two roots of the quadratic 

indicial equation + (6 - a - l)s + 3/4 - 6/2 = =^ s± = ^(-^5 ± a/21). Contrary to 
what happened above for the integration constants entering the homogeneous solutions in 
Ai and i^i, the /^i-integration constants cannot a priori be physically ignored. However, if 
we follow the spirit of [27] and of the leading requirement (64), i.e. if we try to construct the 
simplest vacuum solution which is entirely determined by the GR-like integration constant 
Rs, we can discard the above homogeneous solution and continue the iteration by selecting 
for fii the particular solution determined by the inhomogeneous "source terms" in the 
/ii-equation. This solution has the following solution: 



Note that the exponents s+, s_ entering the //i-homogeneous solution are irrational. Taking 
into account the possibility of adding this homogeneous solution would introduce a new 
sequence of characteristic length scales (25) which would mix in a complicated way with 
the other (rational) scales. 

Continuing in the same way (always discarding the homogeneous solutions entering the 
II equations) we end up with a double series for the three functions f — z/. A, // of the form: 





(75) 




(76) 



^ n=0 k=0 



'rT^ 



n 




k 



f 



nk 



( 



) 



(77) 



22 



where a = 1 for A, and a = 1/2 for For notational simplicity, we have not indicated in 
Eq. (77) the presence of the logarithmic terms log{R/S) which appear each time the "source 
terms" for the perturbations in A and u contain a power of R which match (modulo a factor 
R"^) the power of i? of a homogeneous solution of the left-hand-side. The arbitrary scale S 
that one can introduce in these logs is not physically important because a rescaling of S 
corresponds to the addition of an (ignorable) homogeneous solution in A and u. [Because of 
the irrational character of s+, s_, no logs primitively enter from the perturbation equations 
for the fj, variable] Without loss of generality, we can then choose S = Rs as scale in all 
the logs. 

The final result (77) thereby involves only two length scales Rs and the length scale 
entering the monomial ^^^1/2 ■ For this expansion to make sense, we should make sure that 

both ^ and "^^^^^^ be (much) smaller than 1. This requirement is satisfied if we look at 

distances Rs <^ -R <^ Ry, where Ry is the length scale introduced in Eq. (3), which, using 
the general definition (25) can be written as: 

Rv — Ra,\ (78) 

To summarize so far: we have verified part of the claim of [27], namely the existence of 
a (particular) double expansion, a priori valid in the intermediate range Rs <^ R <^ Ry, 
that (i) represents a formal (vacuum) local solution of the Einstein equations modified by 

the mass term (65), (ii) is continuous as 0, and (iii) satisfies Eq. (64). 

The weak points of this argument are several. First, the specific local solution chosen 
in Eq. (77) corresponds to a postulated energy momentum tensor which does not stem 
from an action, and we are going to see that a similar solution does not exist for the action 
used as a starting point in [28]. But the weakest point concerning the expansion (77) is the 
lack of rationale indicating that this local solution does approximate the unique global one 
satisfying the standard boundary conditions at infinity (16) and being regular at the origin. 
The detailed argument above, showing that two arbitrary integration constants have been 
set to zero to obtain Eq. (77), does not make it a priori probable that Eq. (77) happens to 
satisfy all the required boundary conditions. 

Before tackhng the more complicated "global" issues, let us first study the issue of the 
existence of formal Vainshtein-type expansions for action-based Held equations. We start by 
considering the potential V*^^) that appeared in the recent paper of [28] , which re-discussed 
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the result of [27]. Trying to repeat the same procedure as above we face the following 
problem. The equation for /iq now reads: 

Y^~2R~^ 4 2~ + '"-° ^^^^ 

The //-quadratic terms in this equation are different from those in Eq. (69). This difference 
is crucial because if we look again for a solution satisfying Eq. (64), i.e. admitting an 



expansion of the form //q = jiQo^J Rg / R-\- Hqi{Rs / R) H we find that there exists no such 

solution^. [Formally, the only such solution starts with an imaginary coefficient: //qo = 



The same problem persists if we consider the potential V^^^. Then the same equation 
reads: 

Ao i^o 2//0//0 /^o /^O/^o , n (Q(\\ 

:r^"2^"~^"T"^ + --- = ° 

and would again formally give an imaginary solution: //qo = '^'^\f\- 

The situation is better, however, for the potentials V^^^ and V^^^. The approximated 
equation for V^^^ is then: 

Ao T^'o , 2/^0/^0 , /^o , /^o/^o , n /'qi ^ 

^"iR + ~^ + T + ^ + -"-° ^^^^ 



which gives //qo = -y^ and for V^^^ we have: 



Ao QfJ^ofJ'o I 3/^0 3/io/^o , n rQo\ 

:r^"2^ + ~^ + ^ + ^ + -"-° ^^^^ 

which gives //oo = \ r^- Thus, at least in these Vainshtein-type expansion similar 

to (77) exists. Again, for these expansions to make sense one should be at the intermediate 
range of distances Rs -C -R <^ -Ry, with Ry given by Eq. (3). 

Let us note in passing that the new scale selected by a Vainshtein-type construction is 
not at all universal for different mass terms. For example, let us consider the potential [11]: 

V^''^ = {f9f'\<y2-<yl) (83) 



^This contradicts the statements of [28], which were supposed to concern the mass term V^^^. They, 
instead, seem to correspond to the postulated T^lf ^ of [27] . 
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where (T„ = Ea(lnAa)"', with Aq denoting the eigenvalues of f -"^g (see below). Then, the 
non-linearity expansion of the constraint (35) gives: 

If we introduce the short-hand notations e = mRg and x = R/Rs, we end up with 
instead of Eq. (77), expansions of the type: 



n=0 k=0 

with a = 1 for A, and a — 1/3 ior /i. This expansion has the same form even if we add the 
stabilizing A(t| term (see [14]). By the same argument as before, this gives the following 
range of distances for which the above expansion makes sense: 

A final example would be to consider the following potential: 

V(°) = ifgf'al (87) 

for which the mass of the graviton is zero, but still we have a non-trivial potential which 
breaks general covariance. Then the expansion at lowest order looks like: 



f = x ^ [foo + fioeHogx + ---] 

for all three functions /i, v and A. The formal range of validity of this expansion is now: 

Rs^-R-^ Rse^'^"^^^^^ (89) 

At this point, it is important to remark that, contrary to what one might think (in view 
of the controversy over the "discontinuity" issue) there exist many local solutions of the 
massive-gravity field equations which are continuous as — > and which, therefore, might 
suggest that massive gravity is compatible with local tests of GR. The specific Vainshtein 
construction is (when it works) one way of exhibiting such solutions (constrained to the 
obtention of strictly stationary solutions). However, we have seen above that this specific 
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way fails in important cases, such as the simplest mass term V^^^. However, even when 
there does not exist a /iq satisfying the simple requirement (64), there might still exist 
other solutions. Indeed, the equation that /lo has to satisfy is a second-order quasi-linear 
ODE and one can always construct local solutions (depending on two arbitrary parameters) 
for such an equation. [The problem then is to show that the domain of existence of the 
solution is large enough to cover the phenomenological tests.] In fact, if we forget for a 
moment about the restriction to strictly stationary solutions, Ref. [11] has indicated how 
to construct general classes of solutions, evolving on the time scale m^^ (thought of as the 
Hubble time scale), which are O(m^) close to any desired GR solution everywhere in space. 
All these constructions leave open, however, two basic issues: (i) can these local solutions 
be extended to globally asymptotically flat solutions, or, if they cannot, (ii) can they be 
extended to a cosmological solution which is a natural "attractor" of the cosmological 
dynamics, so that it is indeed natural to use such a solution to describe the universe around 
us? 

Before tackling the issue of whether the formal, particular Vainshtein-type expansions 
are part of a globally regular solution (both at the origin and at infinity), we study in the 
next section the structure of the field equations in a different gauge. 

6 The c, c, b variables 

Let us consider the a,b,c gauge (27), and replace the three basic variables a,b,c by the 
equivalent combinations c, c, b where we define: 

b^ab , c^ca^ (90) 

Let us indicate the link between the X, /i, i/ variables and the a,b,c ones (or the equiv- 
alent c, c, b ones). We first compute from the three functions a(r), b{r), c(r) the functions 
A(r), /x(r), u{r) as follows: 

a2(r) = e^W , ^^^r) = e^^+^M (^1 + , c2(r) = rV^ (91) 

Then we make the coordinate transformation r — > by inverting R = c(r) = re^^"^^^^. 
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The advantage of introducing these variables is that they simphfy the Lagrangian for- 
mulation of the dynamics. Indeed, the action can easily be checked to be: 



where we have introduced V = V/ sin^^. Note that V has the scaling behaviour V = 
r'^f{c/r,c/r,b). We have denoted ' = d/dr. The pair (c, c) appears here as "light-cone" 
coordinates for a "relativistic particle". We refer to the cosmological studies of massive 
gravity [11,14] for a discussion of similar "relativistic particle Lagrangians" . Note that the 
variable 5 is a "gauge-like" variable {radial lapse) which does not have a kinetic term, just 
as the (relative) time lapse variable e^ did in cosmological studies. Therefore the equation 
of motion for b is an algebraic equation, instead of the second-order ODE's one gets for the 
dynamical variables (c, c) . 

The equations of motion for this gauge choice are: 

bdr\bdrj 8 b dc c^*^ ^' ^ ' 

1 dcdc , dV „ ^ 2t-. /^^n 

A consequence of these equations are the conservation equation of the matter energy- 
momentum tensor: 

dr dr dr 

Note that the structure of the equations of motion in vacuum {p = Pr = Pt = 0) 
is relatively simple: Eq. (95) is an algebraic equation for b whose solution yields b — 
B{dc/dr, dc/dr, c, c, r). Inserting this expression in the first two equations then yields two 
second-order, non-linear ODE's for the radial evolution of c, c. 

We note that it is sometimes convenient to reformulate the field equations by changing 
the radial variable from r to the radial analog of the "proper time" , namely f defined by 
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d-r = bdr. Then the set of field equations read: 

14 « 

Note the appearance of a fourth equation in the system hnked to the definition of f . The 
pressure-balance equation takes the same form as Eq. (96) with the replacement dr — > dr. 

The use of the variable f corresponds to writing the two metrics as: 

ds^ = -a\f)df + + c\r)dn^ (101) 
a''{r) 

dsl = -dt^ + + r^{r)dn'' (102) 



6.1 The General Relativity m = limit 

To get some famiharity with the c, c, b variables let us see how one can derive the 
Schwarzschild solution in this new language. In GR 6 is a gauge variable, and we can 
set 6 = 1 for the exterior solution which implies that, after fixing an integration constant, 
we can write r = r. Alternatively, we can start from the f-formulation and set to zero the 
terms 0{irP'). This yields the very simple evolution equations d'^c/df'^ = = d^c/df^, sub- 
mitted to the constraint {dc/df){dc/df ) = 1. The general solution is c = ar + f3,c = 'jr + S, 
with the constraint a'~f = 1. The solution depends on three arbitrary constants. However, 
we can use the remaining rigid gauge symmetries of the problem (f' = af + b,t' = ct) to 
set q; to 1 and f3 to zero. This leaves only one physical integration constant, 6 = —Rs, in 
terms of which the solution finally reads: 

c — f — r , c — f — Rs — r — Rs (103) 



For the interior solution, and for a constant density star, we choose the gauge variable 
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to have the form: 



1 - 2 

so that we continue having c — r. Then the c function reads: 



1 - ^ 1 

^0 (104) 



2 



3 I Rs 1 /, R 



Note that the interior solutions can be obtained from the ones of the A, /j,, u gauge [Eqs. 
(39), (41)], using the connecting relations (91) with /j, — 0. 

6.2 The perturbative limit for m 

We can look again how the results change in the linearized approximation for the massive 
case in terms of the c, c, b. In this case we linearize as: 

c^r + 5c , c = r + 5c , b = 1 + 5b (106) 

We will only present here the exterior star solution to verify the correspondence with the 
A, u gauge. The vacuum equations (97), (98), (99), (100), using the r variable become: 

2 

5c" - 5b' = ■^{r5b + 5c - 5c) (107) 

2 

_ Tfl — 

5c" - 5b' = —{Sr5b + 35c + 5c) (108) 

9 

_ 777 

5(! + 5^ - 25b = —r{35c + 5c) (109) 

for all four initial potentials (5), (6), (7), (8). From Eq. (109) we can solve for 8b and 
substitute back to Eqs. (107), (108). Adding the resulting equations we obtain: 

r{5d - 5t) + (1 + mV^)5c - (1 - 3mV^)(5c = (110) 

Differentiating the above equation and substituting the difference 5d' — 5d' from Eq. (107) 
(with 5b substituted), we obtain a second independent first order equation: 

r(25c' + 5d) + ^(2 - mV)5c+ ^(22 - 3m\'^)5c = (111) 
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Prom Eqs. (110), (111) we can solve for 5c and substitute to any of them. In this way we 
get a second order differential equation for 5c: 



5d' + .^^ , ,. fc' - ^±^5c = (112) 
r(2 + m^r^) H 

while 5c and 5h are then given as: 

At — 1 777^ 

5c = -35c - —5(^ , 5b = -5c' + 5c' - -pr(35c + 5c) (113) 

Thus, although the system (107), (108), (109) at first sight seems to involve four inte- 
gration constants, in fact only two of them are independent. Finally, we can integrate Eq. 
(112) and write the general solution decaying at infinity as: 



5c = ^ fl + — + ) e— ^ (114) 

4 \ mr (mrY 



5c = ^ f-3 + — + e-'- (115) 

4 \ mr (mry J 

5b^-—(l + — + e-"^" (116) 

2r \ mr [mry J 

where we have set the second integration constant corresponding to the exponentially grow- 
ing solution to zero. It is easy to see, from the connecting relations (91) and after defining 

R = re^^'^, that these solutions are identical with the ones obtained in the A, /i, v gauge 
[Eqs. (47), (48), (49)], with the same integration constant Ci. The latter, as we saw in 
section 4.2, for Rq -C m~^ is Ci ~ ^Rs- 



7 Singularity in the general asymptotically flat solu- 
tion 

Let us now discuss the central question of whether non-linear effects do cure or not the 
"discontinuity" (as — > 0) with respect to GR exhibited by the linearized approximation. 
We recall that the specific conjecture made by Vainshtein in this respect was that the 
non-linear "dressing" of the unique perturbative solution (defined by resumming Poincare- 
covariant diagrams, as sketched in [28]) leads to a solution of the type constructed in section 
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5. There are several ways one can take to investigate this conjecture. Either one can start 
from the center of the star with the requirement that it is locally flat and integrate outwards, 
or start from "infinity" (practically from some distance ^ m."^) with the requirement of 
asymptotic flatness and integrate inwards. In addition, we can work in several different 
gauges. 

Our work had to be necessarily numerical because it is a difficult dynamical question to 
control the global structure of a non-linear theory such as massive gravity. Though we have 
done simulations both ways (inwards and outwards), and with several gauges (A, /x, u, and 
a, b, c), we shall mainly report on the results obtained by integrating inwards, in the a, b, c 
gauge. This way of proceeding has indeed two advantages: (i) the solution is essentially 
unique (once the equation of state of the matter is fixed) in the sense that its "initial state" 
at infinity is physically determined by giving oneself only one dimensionless parameter, 
mRs, and (ii) the (inwards) radial evolution of the convenient variables c, c, b is similar to 
the cosmological evolution studied in [14], which allows us to draw some intuition about its 
qualitative behaviour. Note that, by contrast, starting the evolution from the center obliges 
one to consider a two parameter family of solutions. One then needs to "shoot" from the 
center with a one-parameter family of initial data until one eventually matches the unique 
exponentially decaying solution at infinity. 

Let us start by discussing the exterior region. [As we shall see this suffices to conclude 
negatively.] The dynamics of the c, c, b variables is given by the Lagrangian: 

We start at infinity with the solution (114), (115), (116). For this solution, the "kinetic 
term" c c is positive (and near + 1), which means that the analogous "relativistic particle" 
initially moves on a spacelike worldline. Let us recall that, as b has no kinetic term, the value 
of b at any "moment" {i.e. at any radius r) is determined by extremizing the Lagrangian 
(117) considered as a function of b. As long as the sign c c > does not change, the first two 
terms in the Lagrangian (117) would suffice to "confine" b around a unique extremum ~ 

Vd~d. If ever the dynamics tends to decrease the value of c c, this will tend to drive b towards 
smaller values. In fact the dynamics is due to the interplay between the "kinetic terms" 
in (117) and the potential terms and we must also take into account the 6— dependence of 
V (b, c, c, r) . This dependence is somewhat complicated, and has a form which depends on 
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the specific mass term considered. For our four fiducial mass terms V^*\ i = 1,2,3, 4 we 
have: 

V(i)=Pf(i) + 6G'(i) (118) 

V(3) = ^F(3) + 6G(3) (120) 
b 

■pW = l^F^^^ + G(^) (121) 
b^ 

where, F^^^ = r'^F^'^^ (^, ^), G^*^ = r'^G^^^ ^). However, this structure is not enough for 
concluding about the "confining properties" of the potential terms, because the coefficients 
of the various powers of b above do not have a definite sign. Indeed, the explicit form of 
these functions is: 

= ^F^^) = (-12 + 12^ -2-^+6^- 4^] (122) 



= = cM -2 + 6- - 4^ (123) 
J.I \ c r^c J 

= = ( -12 + 12^ - 2^ + 6- - 4- (124) 
r'' \ c cc ) 



We can, however, have an analytical idea of the natural tendency of the dynamics by 
looking at the linearized solution (114), (115), (116) which results from the combined effect 
of kinetic and potential terms. If we look at Eq. (116) (remembering that C\ > 0) wc see 
that all the terms in the variation of b arc negative, i.e. that b tends to decrease from its 
initial value 6=1. In addition one can deduce from Eqs. (114), (115) that the crucial 
quantity c c, after increasing, starts decreasing when r^m~^. As we said above, such a 
decrease of c c tends to further drive b down. 

If we first consider the mass terms i — 1,2, we see that they contain only positive 
powers of b. Whatever be the sign of the coefficients of these terms, such functions are 
quite inefficient in preventing b from decreasing all the way towards zero, if the tendency of 
the "confining" term c c/6 in (117) is to drive b there because c c happens to decrease below 
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zero. Therefore, in the cases i = 1, 2 it might a priori happen that the dynamics drives b 
toward zero at some finite radius, similarly to what we found in the cosmological study [14] 
where mass terms with weakly confining features tended to drive the gauge function to 
zero in a finite time. The numerical simulations we performed for the potentials i — 1,2, 
with various values of mRs (in the range 10~^ < mRs < 10~^) have shown that indeed 
such a behaviour is generic: we systematically find that the inward radial evolution ends 
up in a singularity^ where the variable b tends to zero, while at the same time c c ^ 0. 

We have also run similar simulations with the other mass terms i = 3,4, and we found 
the same singular behaviour: 6 ^ at a finite radius, while, at the same time, c c ^ 0. 
In the latter cases, the apparent "potential barriers" F'^^^b and F^^yW' that might have 
helped to prevent b to tend to zero turn out to be ineffective, because the (non-positive 
definite) coefficients are found to be driven by the dynamics towards zero at the same time 
that c c — > 0. 

Our final conclusion (based on our numerical simulations) is therefore that for all values 
10~^ < mRs < 10~^ (and probably also for all smaller values mRs < 10~^) the unique, 
asymptotically fiat SSS solution (which is well defined near infinity) develops a singularity 
at some finite radius. The latter radius depends on the mass potential and is roughly 
between Rp and m~^. 

To complement our numerical study, we also performed simulations (still in the a, b, c 

gauge) that start at the center of a star. This case is more delicate because generic (regular) 

initial data at r = depend now on two arbitrary parameters, and we expect that only 

a "line" of data can evolve into an asymptotically flat solution. If we consider a weakly 

self-gravitating star, it is relatively easy to choose initial data at the center (with b, c and c 

all near one) that satisfy the constraint (95). Then, when considering the evolution system 

with respect to the "proper radius" f, general theorems on the continuity (with respect to 

small parameters appearing in the coefficients) of solutions of ODE systems guarantee that, 

if m^i?0 is small enough, any regular interior solution in the GR limit will be smoothly 

deformed into some massive gravity interior one. This provides a way to exhibit classes 

of local solutions of massive gravity which are continuously connected to GR solutions. 

However, a hard dynamical question (which is not covered by the local continuity theorems 

for ODE's) is to know on which length scale these solutions stay close to GR solutions, 

'^Within the context of massive gravity the value 6 = corresponds to an invariant singularity of the 
bi-gravity configuration (g,f) because one of the eigenvalues of f~^g is 1/6. 
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and what are their asymptotic behaviour at infinity. We have addressed this dynamical 
question by numerical simulations (in the a,b,c gauge). We found that all our simulations 
evolved a singularity at a finite radius (where b run away toward zero). This singularity 
occurred at distances between Rp and depending on the mass term. In view of this 
singularity, we could not address the issue of "shooting" the initial data with the aim of 
matching an exponentially decaying solution at infinity. 

Wc have also performed simulations in the A, fi, v gauge. Their results have again all 
be negative: all solutions starting from the center developed some strange behaviour for 
large radii, and never showed any tendency to match the perturbative solution at infinity. 
Let us describe some of our results in the A, /x, v gauge. 

In view of the nature of the ODE system in this gauge we a priori need to give initial 
data, at 0, for A(0), ^^(0), /x(0) and /i'(0). First of all, we know that in order not to have 
a conical singularity we should have A(0) = 0. Second, if we remember that the function 
// defines the hnk between the "fiat" and the "curved" radial coordinates r = it!e~''(^^/^ , 
the regularity (in Cartesian-like coordinates) of this link implies that we must also require 
//'(O) = 0. Finally, we have only two initial data: say C = i/(0) and D = //(O). 

The leading terms for the three functions are: 

A = AR^ + ... (126) 
1/ = C + BR'^ + ... (127) 
D + ER^ + . . . (128) 

In GR the quantities A and B are related to the density po (for a constant density star) 
and the central pressure Pc as: 

?>Agr = SttGpo (129) 

2Bgr - Agr = StiGPc (130) 

while the constant C is arbitrary. In the massive gravity case, we obtain the following 
relations: 

3A = SnGpo + m^Ei (C, D) (131) 
2B - A = 87rGPc + m^j:2{C,D) (132) 

for all four potentials V^*), i = 1,2,3,4. The functions Si(C,D) and ^2{C,D) are zero 
when we sit on the point C — D — 0, but they may have non-trivial roots (we will see more 
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in the next section). In order to study the conjecture of Vainshtein (which assumes that 
we are near GR around the source) we explored the region of the (C, D) parameter space 
that is in the neighbourhood of (0, 0), namely -0.1 < C < 0.1 and -0.1 < L> < 0.1. 

For such a choice of initial data, one finds that, for small enough m, the solution stays 
very close to the GR result until a radius which is very roughly of order Ry and then the 
solution starts to deviate considerably. Note that these numerical simulations provide still 
another way of constructing local solutions of massive gravity that are close to GR ones. 
Note also that, when using a matter density model where p undergoes a finite jump at 
R = Rq, the variables A', u" and n" suffer corresponding, correlated jumps. 

However, even when changing in a continuous manner the initial data, we never found 
any solution which, at large radii, tend to match the pcrturbative one. In fact, all solutions 
tend to run into a numerically unstable behaviour, where the simulation stops after a finite 
radius. We saw no tendency that the functions at distances much bigger than Rp match 
the pcrturbative solution ones. We could not run the simulation far enough and could not 
understand if this denotes that there is a singularity at some finite distance from the star. 
However, by numerically transforming the A, fi, v variables into the c, c, h ones we have 
found that at the point of numerical difficulties of the A, /U, v gauge, there was a tendency 
for the c, c, h variables to exhibit their usual singular behaviour (6^0, c c ^ 0). 

Summarizing: Our numerical simulations have falsified the conjecture put forward in 
[27] and [28]. Far from helping to regularize the discontinuities and attendant C(m~") 
terms exhibited by the linearized approximation, the non-linear effects in massive gravity 
aggravate the situation by developing singularities at a finite radius when starting from the 
unique, asymptotically fiat solution. 

This negative result closes, in our opinion, one possibility. However, it still leaves open 
several possibilities for eventually reconciling massive gravity with phenomenology. Indeed, 
our study so far has focused on only one type of possible asymptotic behaviour at infinity: 
namely, (16), i.e. the case where the physical metric asymptotes the "trivial" vacuum 
solution. In the next section, we shall point out that even the simple-minded mass terms 
(5)-(8) might admit other (translationally invariant) vacuum solutions. Finally, in the last 
section we shall consider the possibility that continuity be restored for solutions which 
asymptote some non-trivial cosmology at infinity. 
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8 Non-trivial translationally-invariant, spherically sym- 
metric vacua 



Returning to our original potentials, let us examine all their possible translationally in- 
variant, spherically symmetric vacua. A translationally invariant vacuum is a solution (in 
absence of matter) for which, in some gauge, both g^^^, and /^i^ are constant {i.e. indepen- 
dent of the spacetime coordinates x'^). Requiring spherical symmetry will impose a further 
constraint on the symmetry group fixing the bi-metric configuration {g^u, f^u)- 

The first constraint translates in the statement that the matrix g^ = f^^gua can be 
globally diagonalized as diag(Ao, Ai, A2, A3), where the eigenvalues Aj are constant. In other 
words, there exist coordinate systems where the two metrics read: 

f = -dt'^ + dx^ + dy^ + dz^ , g = -Xodt^ + Xidx^ + Asdy^ + X^dz^ (133) 

Requiring spherical symmetry then imposes the constraint Ai = A2 = A3. This leaves only 
two constants. 

Finally, rewording the result in terms of the X, /i, u gauge, one easily finds that the 
above requirements are equivalent to the constraints: 

A = , e'^ = Ai = const. , = Aq = const.' (134) 

Technically, we are then left with looking for extrema of the mass terms with constant 
values of A, /x, u of the type (134). 

For the last three mass terms we found that they only admit the trivial vacuum (Ai, Aq) = 
(e^e'^) = (l,l). 

On the other hand, we found that the first mass term V*^^^ admits three possible vacua. 
This can be seen by looking at the Xli(C, D), Xl2(C, D) functions in Eqs. (131), (132): 

Ei(C, D) = -^(e^ - l)(e^ + 3e^ - 2) (135) 

E2(C, D) = i(5e^(e^ - 3) + Te^^ - 3e^ + 6) (136) 

remembering that v — C, /i — D in vacuum. Then the conditions Ei(C, D) = 0, 
T,2{C,D) — have three solutions: the trivial vacuum (Ai,Ao) = (e'*, e'') = (1,1), plus 
two new ones: 

(Ai, Ao) = (e^ en = , (Ai, Ao) = {e^, e^ = Q, (137) 
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Note that the first non-trivial vacuum corresponds to g^i, — jf^u- Such a conformal situa- 
tion respects not only spherical symmetry, but in fact also the full Poincare symmetry. By 
contrast, the last vacuum breaks the Poincare symmetry to a mere Euclidean one. 

Let us examine now the perturbations around these vacua for the first potential. For 
the trivial vacuum, we already know that the linearized perturbations feature a Pauli-Fierz 
mass term, i.e. a ratio a = 1 between the two terms, with a spin-2 mass, say 1712, identical 
to the mass parameter entering the action, i.e. m\ = m^. For the (|, i) vacuum, one finds 
that a = — I and m\ = 16m^. The fact that we do not have a Pauli-Fierz mass term 
means, as we recalled above, that we have a ghost-like scalar field (h). In addition, the 
fact that the value of a is not between 1/4 and 1 means that this ghost degree of freedom 
is also tachyonic, since itIq — — |m|. As we said above, such excitations are pathological, 
and we do not wish to consider them (though we know that this is a cheap way of ensuring 
continuity as — > 0). 

Finally, the third vacuum (|, |), does not give an equation of the type (17), but rather 
instead of the h^i, — arj^yh mass term, one finds the following contribution: 

D;^Kx + E^,h (138) 

where the tensors D^'^ and E'^i^ are constructed from the two background metrics fjf^i, and 

The important point is that these tensors for the background we are considering are not 
proportional to 5^5^ + ^ v) and f]^j^y respectively. This happens because the "metric 
condensate" is not proportional to the background of the observable metric, and therefore, 
as we said, breaks the underlying Poincare invariance. Thus, as we see from Eq. (138), the 
equation of motion for the perturbation h^^ will have different structure from Eq. (17). We 
leave to future work the task of redoing, in this new setting, the analysis of [17]. 

In the next section, we shall finally study other types of vacua, which admit a cosmo- 
logical interpretation. 
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9 Continuous "non-co-diagonal" solutions with de Sit- 
ter asymptotics in massive gravity 



In this section we explore another possibihty for eventually reconciling massive gravity with 
phenomenology. This possibility is based on an interesting type of spherically symmetric 
stationary solutions of massive gravity found long ago by Salam and Strathdee [39] . [These 
types of solutions were then generahzed to a full bigravity setting [40].] These solutions 
tend, at spatial infinity, towards a de Sitter solution, instead of one of the translationally- 
invariant vacua considered above. 

The aim of this section is two-fold: (i) to describe the "universality class" [11] of mass 
terms which admit such solutions, and (ii) to show how the vacuum solution of [39] can 
be extended to a globally regular spacetime representing a "star" embedded in a de Sitter 
background. 

9.1 Spontaneous symmetry breaking in massive gravity 

Up to this section, we have only considered spherically symmetric solutions which could 
be represented in a co-diagonal way, i.e. such that both g and f are diagonal in the 
same coordinate system. As we indicated, perturbation theory (around the trivial vacuum) 
necessarily select such a type of solution, if the matter sources are supposed to be "at 
rest" . However, non-perturbative solutions might introduce an interesting twist in the link 
between the two metrics g and f . A convenient formalism for describing this possibility 
consists in introducing a "hnk field" Y{x) (also called "Stueckelberg field") relating the 
two metrics g and f. This formalism first appeared in [44] and, in more detail in [45], 
and has been recently developed in [46]. More precisely, instead of writing the theory of 
massive gravity in a coordinate system which is common to both metrics, and of writing 
the mass term in terms of the difference g^u{x) — f^u{x) between the two metric tensors 
(supposed to live on the same manifold, and expressed in the same coordinate system), one 
can explicitly introduce the map Y between the g-Riemannian manifold, and the f-one. If 
we use independent coordinate systems on both manifolds (say x^^ on the g manifold Aig, 
and X" on the f-one Aif) the abstract map Y takes the explicit form X" = y"(a;^). This 
map (from Aig towards Ai /) "pulls back" the metric f into an image metric Y^i which lives 
on A4g. This leads to expressing the mass term in terms of: 

H^^Ag, f, y] = - d^Y-d,Y^Up{x) (139) 
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Then the mass term (and therefore the action) becomes a functional of three fields: g, f 
and Y. In the massive gravity context, f is a non-dynamical flat metric, and we can always 
choose to express it in Lorentzian coordinates, i.e. with fap{x) — rjap. Then the action 
depends only on two dynamical fields: g and Y. The usual formalism (used above) consists 
of using the gauge invariance of the theory to fix y = id (which destroys the explicit 
gauge invariance) . Ref . [46] advocates the usefulness of not using this "unitary gauge" , but 
instead of working with an explicitly gauge invariant theory with two fields: g and Y. The 
immediate problem with this proposal is then that the kinetic terms for the Y field are 
very non- linear because they result from replacing (139) in the mass terms. [For instance, 
the existence of a sixth degree of freedom is hidden in the highly non-hnear nature of this 
kinetic term.] 

One can view F as a gravitational analogue of the Higgs field. Ref. [46] has indicated 
how to deal with perturbation theory around the "trivial Higgs configuration" Yf^^x^^) = x", 
by expanding the theory in powers of the "Goldstone" field n°'{x) such that Y°'{x'^) = 
x" + 7r"{x). The analogue, in this language, of the perturbative co-diagonality of two SSS 
metrics is the following. If we consider some static matter source (in M.g), perturbation 
theory for the Goldstone field will never generate a time-component Tr^{x) (which induces 
a non-zero mixed component goi{x) in the physical metric). 

In this language, one needs to be in a non-perturbative situation of "spontaneous sym- 
metry breaking" to expect such a 7r°(x) to develop. This situation was in fact recently 
considered, mutatis mutandis, in [14]. There it was noticed that besides the stable family 
of bigravity cosmological solutions, where the two metrics gi , g2 can be written in a time- 
orthogonal gauge {goii = = ^'012) and where the remaining metric coefficients {gooydij) 
depend only on time, there might also exist, under some conditions, families of "symmetry 
breaking" solutions where the relative "shift vector" 6* (proportional to (72* in the gauge 
where g^n = 0) is not equal to zero. The possible non-zero values of 6* were the extrema 
of the mass term V^(b) ~ Vq + ah^ + feb^ + 0{h^). Similarly to the Higgs mechanism, for 
certain shapes of the function V^(b) {e.g. when a < and 6 > 0) there might exist, besides 
the trivial extremum b = 0, non-trivial extrema (b 7^ 0) of V^(b). The precise condition for 
this to happen was studied in the Appendix B of [14]. Wc are going to recover the same 
condition for the possibility of a similar symmetry breaking in the context of SSS solutions. 

In this context, the co-diagonal form (28), (29) is not the most general form compatible 
with spherical symmetry and stationarity. The most general form depends on 4 functions. 
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and can be written in the following A, a gauge: 



dsl = -{dT + adRf ' e-^^^'^dR^ + e-^^^""^ RH^'' (141) 

Here all the functions A, /x, z/, a depend only on R. The above gauge has the advantage 
of separating the variables into two groups: the "physical" variables A, v (which are di- 
rectly observable in gravitational experiments), and the "gauge" variables /i, a which enter 
only the unobservable background fiat metric. The gauge variables enter the action (4) 
only through the mass term — (m^/4)V(f~^g). Therefore one can immediately obtain the 
equations of motion of o; by varying only the mass term: 

i-V(f-ig) = (142) 

I^V(r^g) = (143) 

Equation (142) above generalizes the constraint (35) to the case where a ^ Q. Let us 
consider here the new equation (143). A crucial point is that a enters V only algebraically. 
Therefore, similarly to what we recalled above concerning the dynamics of the cosmological 
shift vector b, non-trivial SSS solutions will exist only if V{a) admits non-trivial extrema 
with respect to a. To understand what this imphes for the "universality class" of such 
mass terms, let us relate the a dependence on the eigenvalues of the gravitational energy- 
momentum tensor TjfJ . Let us recall (from [11,14]) the following bi- g eometrical ia,cts: The 
bi-gravity configuration (g, f) defines, at each spacetime point, a preferred moving frame 
with respect to which g = diag(— Aq, Ai, A2, A3) and f = diag(— 1, 1, 1, 1). [Note that this 
preferred frame is not proportional to the coordinate frame used in the normal gauges such 
as the A, /x, u, a gauge.] The eigenvalues Aa of f~^g depend on R. They can be easily 
computed from the metric coefficients A, /x, u, a. For instance: 

A2 = A3 = e''(^) (144) 

Let us focus on the a dependence of the eigenvalues. Note that a disappears when calcu- 
lating the determinant of f . Therefore the product: 

A = AoAi (145) 
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which is the ratio of the two 2x2 determinants of the (T, R) sub-blocks of the two metrics, 
does not depend on a. Indeed, 

On the other hand, a enters the trace of the (T, R) sub-block of g~^f . Calculating this 
trace yields: 

l + l = e-" + e-*[(l-^)^-„^l = ^ (147) 

Ao Ai / AqAi 

When varying only a, one therefore finds that A2 = A3 and AqAi do not change, but 
that Ao -l- Ai varies proportionally to aSa. On the other hand, the variation of the mass 
term S'^^^ oc —V in the action is related, by definition, to the logarithmic changes of the 
eigenvalues by 5S^^^ — (l/2)y^T(f)"51og Aq, where T^^^^ denote the eigenvalues (in the 

preferred frame which is easily seen to diagonahze also T^^) of the gravitational energy 
tensor. The final result is that the variational derivative (143) is found to be proportional 
to the product: 

^ ; , ' (148) 

Ao — Ai 

As we said above, a universal solution of this ct-constraint is the trivial case a = 0. The 
spontaneous symmetry breaking (SSB) situation corresponds to the case where it is the 
second factor above which vanishes. Using the formulas given in [11, 14] one can easily 
see that the condition for SSB does not depend on the choice of density prefactor in front 
of V = wV, where w can be either \/—f, ^/—g, {fgY^'^ etc... Finally, we can write the 
condition for SSB in terms of the scalar potential ^(Ao, Ai, A2, A3) = V /w as the existence 
of solutions for the equation: 

,(A.,A„A„A3).^(^-5^) = (149) 

For instance, if we consider the class of potentials V which depend only on ui and (T2, where 
(7„ = Ea(ln Aa)", we recover the condition d^^V — derived in [11, 14]. 

It happens that all the mass terms written in (5)-(8) can satisfy the SSB condition. First, 
we note that they contain only two independent "scalar" potentials V . Moreover, the two 
types of scalar potentials can be mapped into each other by the exchange g f , which is 
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equivalent to the exchange Xa <-> 1/Xa- Therefore, it suffices to consider the potential V^^^ 
defined, say, by: 

V^'^ = I {tr[((g - f)rr] - (tr[(g - f)f-']r} = ^E(A„ - If - (5](A„ - 1))^] (150) 

a a 

The computation of the quantity p defined above then yields: 

^(Ac-l) = d- ^ Ac (151) 

c^O,l c^O,l 

where d denotes the space dimension. In our case, it is d = 3, and we can also use A2 = A3 
to write p — 3 — 2X2. 

Finally, we conclude that for the mass terms V*^^-* and V*^^-*: 

p[V(i)] = p[V(')] = 3 - 2A2 (152) 

while the corresponding results for the other ones (obtained by exchanging Xa ^ 1/Aa) 
read: 

p[V(3)] = p[v(^)] = -3 + 2/A2 (153) 

Therefore the SSB condition (i.e. the equation of motion for a 7^ 0) can be satisfied for 
all four potentials. In the first two cases, it fixes the value of A2 , i.e. of e'^, to be simply 
3/2, while in the last two cases one must have X2 — — 2/3. [The corresponding values in 
space dimension d would be d/{d — 1) and {d — i)/d, respectively.] It is a happy technical 
accident (which simplifies the subsequent derivations) that these values are constant. This 
is, however, not crucial. The essential point is that the SSB condition p{Xa) — admits 
solutions. More general types of SSB potentials might admit solutions of p(Ao, Ai, A2) = 
along some curved hypersurface in A^-space, say A2 = /(Aq, Ai). 

9.2 Non-codiagonal solutions in massive gravity 

Let us now study the global structure of SSB solutions (both in the exterior and the interior) . 
There are two ways of proceeding. Either we continue using a "geometrical" approach based 
on the systematic use of the invariant eigenvalues, or we use a more pedestrian approach 
based on writing down explicitly all the field equations in some gauge. Before using a more 
pedestrian approach let us briefiy indicate the results of the "geometrical" one in the simple 
exterior case. 
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Having obtained the value of A2 = we can now get a simple constraint on the product 
of eigenvalues A = AqAi. Indeed, we found above (as a consequence of the SSB condition) 
that the 2 x 2 (T, R) block of the gravitational energy tensor satisfies (in its diagonahzing 

frame) T^^^^q = T*^^)J. This condition actually means that, in any frame, the 2x2 {T,R) 

block of T^v is proportional the 2x2 unit matrix. When going back to the physical gauge A, 

/X, u, a this means that the equation of state of the gravitational energy satisfy Pr^^ — —p^^h 
Using this information in the usual Einstein equations (in the exterior) (31), (32) gives us 
the simple constraint that the product e'^e^, i.e., modulo a constant, the quantity A = AqAi, 
is constant. At this stage we know all the eigenvalues modulo the knowledge of, say, Xi{R). 
To get the radial variation of Xi{R) it suffices to write down the conservation equation (36) 
for the gravitational part of the energy-momentum tensor. This is simplified by the result 

Pr^^ — —p^^^ and now reads 

9«Pj^) + |(P,(^)-P/^^)=0 (154) 

As the components {—Pr^\ Pr^\ pI:^\ Pj:^^) of the gravitational energy-momentum tensor 
are checked to depend only on the eigenvalues, the above equation gives an evolution 
equation for \i{R). Actually, one can see immediately a particular solution: namely Ai = 
A2. [One then checks that this is the only possible solution.] Indeed, the isotropy condition 

Ai = A2 = A3 guarantees that Pj:^^ = P}^\ Then we get a solution of the above equation 
because all eigenvalues are constant (indeed, we know already that A2 = A3 is constant, and 
that the product AqAi is also constant). 

Finally, we end up (in the exterior region) with 4 constant eigenvalues, three of which 
are determined to be Ai = A2 = A3 = 3/2 , or 2/3 (depending on the choice of mass 
term), and the fourth Aq being an arbitrary constant (with AqAi = A = const.). Though 
the set of eigenvalues introduces a dissymmetry between the (eigen-)time direction and 
the space directions, there is a further simplification which comes (again because of the 

SSB constraint) from the result above T^^^q = T^^^\. This condition, together with the 
equahty of the three "spatial" eigenvalues, shows that, finally T^^^l = T^^^] = T^^^l = T^^f^. 
Therefore the gravitational energy tensor is fully isotropic in spacetime: T^^^^ = 
where the constant A is some function of the eigenvalues. Such a gravitational energy tensor 
is equivalent to a cosmological constant A. Therefore we conclude, without calculations. 
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that the only SSB spherically symmetric stationary solutions of massive gravity for the 
potentials (5)- (8) must be Schwarzschild-(Anti-)de Sitter solutions. This was, indeed, the 
result of [39]. 

Let us now generalize this exterior solution to a global solution, including the interior 
of a star. For concreteness, we shall consider the specific mass term V^^), (8). The mam 
message of this section would not be altered if we used any of the (5), (6) or (7) instead. We 
know already, from the previous reasoning, the general structure of the exterior solution. We 
only need to compute now the deviation from Schwarzschild-de Sitter (SdS) in the interior 
of the star. We could do the explicit calculations of the interior structure by working, as 
above, with the geometrical eigenvalues, and the corresponding gravitational energy tensor. 
We shall however use a more pedestrian approach, writing down the equations of motion 
in the gauge defined by the line elements (140), (141). [In Appendix B we also present the 
results following the notation of the original derivation [39].] For simplicity we will consider 
a star of isotropic pressure P. The (T, T) and {R, R) equations of motion, after setting 
— 2/3, then read: 



A' 
R 



V 647rG V 2 



871^ 



-3(A+i/)/2 



9m2 



UnG V 2 



g-3(A+i.)/2 _^ p 



(155) 



(156) 



The conservation of the matter energy-momentum tensor gives the familiar equation: 



V 



(157) 



while from Eq. (142) (which replaces the previous constraint fg — derived from the 
conservation of the energy-momentum tensor generated by the mass term) one gets an 
algebraic equation for a with solution: 



a' 



'(R) = ^e-" [3R{\' + v') - 2(e^ - l)(3e" - 2)] 



(158) 



Note that o?{K) undergoes a finite jump (downwards) as R crosses the radius of the star 
(from the interior) . 
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In the exterior, one trivially finds by adding Eqs. (155) and (156) that A' = lie^'^")' = 
which implies, as we said earlier, that A = |(e^"'"'') = const. = Aq > 0. Then from Eq. 

(155) for example one finds the usual SdS solution: 



Then from Eq. (158) one can directly find a: 



,2 2 9A„ 



The requirement of positivity of then restricts < Aq < 4/9. Note also that at the 
position of the dS horizon, the gauge function a diverges since oo. However, this 

is a mere coordinate singularity of the bigravity configuration f, g as is expected from 
the fact that the eigenvalues of f~^g are regular there, being equal, in the exterior, to 
(Ao, Ai, A2, A3) = (3A0/2, 2/3, 2/3, 2/3). A coordinate transformation which regularizes 
both f and g is dT dT + a{R)dR. [Note that the black-hole solution of [39] is also regular 
at the event horizon by the same transformation, something not noticed in [39].] 

The exterior massive gravity solution above involves, as separate scales, only Rg and 
Aq^^A^. Note, however, that (when Aq ~ 1) the scale Rp (2) enters indirectly the solution 
as the scale where the two separate contributions to = 3^6^ become comparable to 
each other. 

In the interior of the star, one has to solve the closed system of differential equations 
(155), (156) and (157). Note that the contribution of the massive graviton in the interior 
equations of motion differs from a pure cosmological constant. A crucial feature of this 
system is that the graviton mass appears only with a positive power, and in front of 
a lower-derivative term. General theorems on the continuity (with respect to parameters) 
of solutions of ODE systems then guarantees that, if m?R^ is small enough, any regu- 
lar interior solution in the GR limit will be smoothly deformed into some massive gravity 
interior one (satisfying the regularity condition \{R) oc R^ at the center). The only re- 
quirement then that the solution is acceptable is that in the interior be positive. We 
have numerically verified that this criterion is not difficult to satisfy. 

Finally, let us note that the only effect of the mass of the graviton in this kind of SSB 
solutions is the introduction of a cosmological tail with an associated scale ~ Aq'^^A^. Since 
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Ao is an arbitrary dimensionless integration constant (lying in the interval < Aq < 4/9), 
one obtains the peculiar result that the cosmological scale of the theory is undetermined. A 
priori, however, one should expect that Aq ~ and that its precise value should finally 
be determined when trying to match this SSS solution with a cosmological one. Note also, 
that since the regular SSB solutions that we constructed differ from the GR ones in the 
presence of a cosmological constant A only inside the star, one does not have observable 
consequences of the type claimed in [34, 35] that might test the model experimentally. 

10 Conclusions 

We have explored the issue of continuity in the purely massive gravity theory for several 
mass terms, which all reduce to a Pauh-Fierz one in the hnearized limit. We have discussed 
in detail the claim of [27] that the resummation of non-linear effects make the massless limit 
smooth. We showed that the type of local expansion (considered in some bounded domain 
outside the source) proposed in [27] does not exist for the simplest, action-based, mass 
term considered in the subsequent, more detailed publication [28]. However, it exists for 
other mass terms. We showed the double series nature of the expansion that was implicitly 
performed in [27] and noted that different mass terms give different kinds of expansions. 

Then we discussed the global nature of massive gravity solutions. When considering the 
inward radial evolution of general asymptotically flat solutions, we find that they all end up 
in a singularity at some finite radius well outside the source. When, on the other hand, we 
consider the asymptotics of a two-parameter class of solutions which are close to GR in and 
near the source, we also find that they end up in a singularity at some large radius. These 
results were obtained in two different gauges, one of which allows one to qualitatively relate 
the singular runaway of the solution to a similar runaway found in recent cosmological 
studies of massive gravity [14]. 

These results falsify the continuity claim (as it was formulated) of [27,28], and vindi- 
cate the statements of [15-17] that all asymptotically flat solutions of massive gravity are 
discontinuous as m — > 0, and phenomenologically excluded. Let us note, for completeness, 
that to fully vindicate this negative result one should still explore the properties of the new 
types of (asymptotically flat) vacua that we have discovered in section 8 above. 

In addition, one should also explore the effect of using "more conflning" potentials, such 
as the potential V oc a2 — o^' -|- Xa-i considered in [14]. There it was shown that such a 
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mass potential can cure the runaway of the gauge variable (when A > 0). One could 
hope that, within our context, it might similarly cure the h runaway by confining h within 
a finite interval separated from zero. However, we do not think that such a modification 
of the mass term is sufficient for curing the discontinuity problem. First, we anticipate 
that any modification of the mass term by additional "confining" terms containing only 
parameters of order unity (such as A in the example above) will have great difficulty in 
confining h because it will have to fight against the extremely large dimensionless num- 
ber h ~ GT/rn? ~ U/{mR)^ which, according to perturbation theory, will tend to enter 
the evolution of the gauge variables b /j ^ h. [And, indeed, our preliminary numeri- 
cal simulations found a runaway of b toward zero with potentials V (x a2 — (Ji^ + Ao"2^.] 
Second, even if one "confines by brute force" the evolution of b by adding sufficiently effi- 
cient terms in we see no reason why the resulting global solution should get close to a 
GR one. We anticipate that, even if the runaway singularity is cured, the resulting (non 
spontaneously symmetry breaking) solutions will be very different from GR solutions, and 
therefore phenomenologically unacceptable. 

Our negative result still leaves open the possibility, proposed in [11,14], that continuity 
and phenomenological compatibility be restored by considering solutions that asymptote 
some non-trivial cosmological spacetime of curvature TZ — 0{'nn?). [Note that this proposal 
differs from what is suggested by the continuity results in (A)dS backgrounds [21-24] which 
need TZ » 0{rn^) to restore continuity.] We have provided evidence for this possibility by 
exhibiting a particular class of "symmetry breaking" solutions with matter sources, which 
generalize black- hole- type solutions found long ago in [39]. As m — > 0, these solutions 
smoothly tend to an asymptotically flat general relativistic spherical star model. The 
"price" for continuity is, however, twofold: (i) asymptotically the solution tends, when 
m 7^ 0, to a de Sitter solution, and (ii) the mass term must belong to one of the special 
"universality classes" dehneated in [11]. The "price" (i) might actually be seen to be a 
virtue, as current cosmological data favor such an asymptotic state for our universe. On 
the other hand, the requirement (ii) is not innocent because it was shown in the latter 
reference that the mass terms coming from brane models do not belong to the needed 
"symmetry breaking" universality class. 

Let us finally note that our positive results (on the rc-cstablishment of continuity for 
symmetry-breaking asymptotically dS solutions) are still far from a satisfactory proof of the 
physical admissibility of massive gravity theories. Many subtle issues need to be checked. 
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such as: (i) can the simple spherically symmetric solutions explored here be deformed into 
inhomogeneous spacetimes able to describe the solar system and the observed cosmological 
universe ? (ii) are such solutions attractors of the cosmological evolution ?, (iii) are they 
stable under quantum fluctuations^^ ? In particular, we worry about the (classical and 
quantum) effect of the "sixth degree of freedom" present in massive gravity [17] whose role 
in the simple solutions explored so far might have been unduly downplayed. 

Let us finally mention that a better setting for investigating these questions (espe- 
cially the "cosmological attractor" issue) would be a full bi- gravity theoij [11,41], which, 
contrary to the purely massive gravity explored here, does not contain anya priori fixed, 
non-dynamical background. In that case, there exist bi-de Sitter "locked" solutions [14] 
which are cosmological attractors (and which are more general than the above considered 
SSB ones). They can be trivially generalized to bi-Schwarzschild-deSitter solutions with 
correlated source masses. An interesting question then is whether these kinds of solutions 
are "spacelike" attractors as well, if we detune the special relation between the source 
masses. 
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Appendix A: The functions ft, Jr and fg 

In this Appendix we will present the functions /(, and fg for the various massive 
gravity theories that we considered in the text. First, these quantities for the postulated 
^°See [47] for recent studies of massive spin fields in (A)dS backgrounds. 
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energy- momentum (see (65)) of [27] read: 



(3 - 26^^ - e") 



(A.l) 
(A.2) 



2R 



1 - e-^ 1 



(l-^)" e^-'^ + (e^ + 2e''-3)(l-^)' 



4 



g-A-M+z^^gi^ + 26^^ - 3) ( 1 



Rli' 



- 6^+^^ 1 - 



-A— /X 



2 y V 2 

i?/x''(e'^ - 3) - 2A' (^1 - {e" + 26^^ - 3) 

-/i' (4(e^-3) + e'^(4 + i?z/')) 

+ 6i?//" - 4e'^i?At" + 26^^(1/' - Rfj,") 



-2e^ + 3 



Rn' 



(A.3) 



For the theories which we have been studying, and for which the equations of motion 
stem from an action principle we obtain the following expressions. 

For the potential V'-^-' (see (5)) we have: 



e^+^(3e" + 2e^ - 3) (^1 - + 3e"(2e^ - 3) + e^(e^ - 6) + 6 



3e^+'^(2e'^ + e'^ - 3) ( 1 - ) + 6^(26^^ - 3) + e'^(e'^ - 6) + 6 



(A.4) 
(A.5) 



f = - — 
8R 



-8(e^ - l)(e^ + e'^ - 3) + 3i?V^(4e^+'* + e'^(e^ + 4) + Ae^ - 3e^ - 12) 



+ ^i?V^(e'' + e"-3)(i?/x' 



6RX'e^{2e^ + e'^ - 3) (^1 



R/j,' 



+ R'^n'u'e^{2e^' + Se" - 3) + 8Rn'{6 + 3e^ - 2e^ - e'^(e^ + 2) - 4e^+^) 



- GRe^+'iu' + Rn") - 2Re\2e'' - + 2,Rii") 



1 - 



(A.6) 



49 



For the potential V^^^ (see (6)) we have: 



f ^ 1 (3i/-3m-A)/2 

jt 2 



,^+1^ 1 



R/i' 



+ (2e^ - 3) 1 



fR = -ie(3^-'^-^)/2(2e'^ + - 3) (l - ' 



(A.7) 
(A.8) 



-(A+3/i+;/)/2 
8 



-2A'e^+''(2e'^ + e'^ - 3) ( 1 - + ^i'^e" + Rv'{2e^ - 3)) (l - ^ 



+ ^(-8(e^ - 3) - AR^i'{e^' + 3) + i?V='(2e'^ + 3) - 2i?V(2e'^ - 3)) 



R 



(8 + 2i?i/' + - - 4) + 2i?V") 



1 - 



(A.9) 



For the potential V^^^ (see (7)) we have: 



e-^ie" + 6^^(36^^ - 2)) ^1 - + 26^(36^^ + 1) - 3e2'^(2e'^ + 1) - e'^ 



(A.IO) 



fR = -^e-^-'^^' 



e^(l + Qe^'ie^ - 1)) - e'^(3e'^ - 2) + e-^(e^(3e'^ - 2) - e'^) ( 1 



(A.ll) 



-A— 2/i 



8i? 



jg/^ -8(e^ - 1) - 2 (^1 - i?A' - 2Rv' + i?Ai'(i?i^' + - 8) - 2i?V 
8(e^ - l)(3e^ - 1) - i?V"(3e^ - 4) + 2i?A'(3e^ - 2) (^1 - 
+ i?//'(3e'' - 2)(8 + - 2i?(3e'' - 2){v' - Rix" 



1-^1 (A.12) 



For the potential V'-"'-' (see (8)) we have: 



ft = -le-(^^+^-+^)/^ (1 - ^ 



e^(3e'^ - 2) - ( 1 



Rfi' 



(A.13) 
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[e-(3e^ - 2) - 



(A.14) 



-(3i/+3A+7;i)/2 

16^ 



2 1 



8 - 6RX' 1 



+ Rn'{3Ru' + 5Rfj,' - 20)^ + 4e^(i?i/'(3e'' - 2) - 4) 



-2 1- 



^ 1^-8 + 24e^ - Rii'iGOe" - 44) 



+ i?V'(15e'' - 14) - 6i?A'(3e'^ - 2) (^1 - 

- 6i?V(3e'' - 2) ) + 16e^(3e^ - 1) 



1_:^^A.15) 



For the potential V^'") (see (83)) we have: 



__p(3^-3m-A)/4 



3/x(/x + i/ + 4) + A(3/x + + 4) 



- 2(2// + 1/ + 4) log 1 - 



Rfi' 



1 - 



Rn' 
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i/(3// + A + 4) + /x(3// + 3A + 8) 
-2(2/x + z/)log fl-^ 



Rn' 
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g-(A+i/+3/i)/4 



32^/l-^ 
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-4(i/ + 2/x) 
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) (3/.' + i.' + A') (^/x (^8 + 3(A + /.) - 41og (l - ^)) 
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4(3Ai'(// + 2v'{v + + 3^^)) 



V ^8 + 3(A + //) -41og (^1 



4i/' 4 + A + 3// - 2 log 1 



(A.18) 



Appendix B: The SSB solution in a different gauge 

In this Appendix we present the SSB solution of section 9 in the notation of the original 
derivation [39] . We now use the following gauge for the two metrics: The background flat 
metric is written as: 

dsl = -de + dr^ + r^dn'^ (B.l) 
while the observable metric is written in the following form: 

ds"^ = -C(r)dt^ + A(r)dr'^ + 2D{r)dtdr + B(r)dif (B.2) 

The off diagonal term D/C is equivalent (after going from the 'fiat" time dt to the "curved" 
one dT — dt — {D/C)dr, and changing dr — {dr / dR)dR ) to the (gauge) variable a above. 

Note the values of the eigenvalues of f in this new notation: 

A = XoXi ^ AC + D^, (B.3) 
Xo + Xi^A + C, (B.4) 
A2 = A3 = B/r^ (B.5) 

For simplicity, we will consider a star with isotropic pressure P. The energy momentum 
tensor is then: 

Tir"^ = (P + PKu, + Pg,.^ (B.6) 
where Ufj, — (— ^,0,0) is the unit velocity along the timelike static Killing vector 
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Cm = (~^' 0' 0)- Explicitly, it reads: 



/ pC 



^^(matt) 



PB 



(B.7) 



PEsin^^ / 



The first step for finding the solution is to impose the SSB condition (149) (vahd when 
a ^ i.e. D ^ 0). For the mass term V^'^\ (8) that we consider for concreteness, this 
yields: 

B 2 



A2 — A3 — 



r2 3 



(B.8) 



Then we can follow the steps indicated above to derive the exterior solution. The exterior 

constancy of A came from considering the effect of the source combination p^^^ + Pr^\ In the 

present notation this corresponds to using AT^^'^^ + CTrr'^^ = 0, to get AT/™'^**-' + CtJ^'^"'* = 
A(p + P), which imphes that ARu + CRrr — 8nGA{p + P). The latter equation explicitly 
reads: 

A'^SttG— r(p + P) (B.9) 
The conservation of the gravitational energy-momentum tensor TjSf' then gives: 

3r A' - 9A2 + (6(A + C) - 4) A = (B.IO) 
while the conservation of matter energy-momentum tensor j'^™'^**) gives: 



P' = 



2~C 



ip + P) 



(B.ll) 



Combining Eqs. (B.9) and (B.IO) we obtain the following expression for the combination 
(B.4): 



A + C^^A+'^-AnG^r'ip + P) 



(B.12) 



One finally needs only one non-trivial field equation to close the system, e.g. the (t,t) 
component of the Einstein equations: 



A2(6 - 327rG'rV) - 4A(C + rC) + ArCA' = Sm^r^VA 



(B.13) 
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The system of radial ODE's (B.9), (B.ll), (B.13) gives, for any given equation of state 
P — f{p), a, closed evolution system for the variables A, P and C. A crucial feature of this 
system is that the graviton mass appears only with a positive power, and in front of 
a lower-derivative term (see (B.13)). General theorems on the continuity (with respect to 
parameters) of solutions of ODE systems then guarantees that, if tti'^Rq is small enough, 
any regular interior solution in the GR limit will be smoothly deformed into some regular 
massive gravity interior one. 

To complete this Appendix, let us give the explicit form, in the A,B,C,D gauge of the 
exterior solution. It is given by A = const. = Aq, A + C = |Ao + |, and by the nontrivial 
equation for C whose general solution is easily found to be: 

C7=^Ao(l-p(r)) (B.14) 



where 

p(r) = ^ + (B.15) 

One then recognizes the Schwarzschild-de Sitter (rather than AdS) solution by changing 
the time variable to: 



df^^ dt-dr^^^^^] (B.16) 

where (3 — (4/9Ao) — 1. The observable metric gets transformed to an exact Schwarzschild- 
de Sitter metric (here represented in a conformally scaled form): 



ds^^l 



[l-p)df^+^dr^ + r^dn^ 
1 — p 



(B.17) 



The connection of these coordinates with the ones used in section 9 is simply R'^ — and 



dT = g^dT. Note also that Rs = J^Rs- 
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